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We report the results of a high-statistics search for H-dibaryon production in inclusive T(1S) and 
T(2S) decays. No indication of an //-dibaryon with mass near the Mh = 2m a threshold is seen in 
either the H —> Apn~ or AA decay channels and 90% confidence level branching-fraction upper limits 
are set that are between one and two orders of magnitude below the measured branching fractions 
for inclusive T(1S) and T(2S) decays to antideuterons. Since T(l, 25") decays produce flavor- SU (3)- 
symmetric final states, these results put stringent constraints on //-dibaryon properties. The results 
are based on analyses of 102 million T(15*) and 158 million T(25) events collected with the Belle 
detector at the KEKB e + e~ collider. 

PACS numbers: 14.20.Pt, 12.39.Ba, 13.85.Rm 



In 1977, Jaffc predicted the existence of a doubly 
strange, six-quark structure (uuddss) with quantum 
numbers I — and J p = + and a mass that is 
~ 80 MeV below the 2mA threshold, which he dubbed 
the /f-dibaryon [1]. An S = —2, baryon-number B = 2 
particle with mass below 2mA would decay via weak in- 
teractions and, thus, be long-lived with a lifetime com- 
parable to that of the A and negligible natural width. 

Jaffe's specific prediction was ruled out by the obser- 
vation of double- A hypernuclei events [2j-|4|, especially 
the famous "Nagara" event that has a relatively unam- 
biguous signature as a aaHc hypernucleus produced via 



S capture in emulsion |3j. The measured AA binding 
energy, B AA = 7.13 ± 0.87 MeV, establishes, with a 90% 
confidence level (CL), a lower limit of Mh > 2223.7 MeV, 
severely narrowing the window for a stable H to the bind- 
ing energy range Bh = 2mA — Mh < 7.9 MeV. 

Although Jaffe's original prediction for Bh — 81 MeV 
has been ruled out, the theoretical case for an i3"-dibaryon 
with a mass near 2mA continues to be strong and has 
been recently strengthened by lattice QCD calculations 
(LQCD) by the NPLQCD 0,0 and HALQCD Q collab- 
orations that both find a bound i7-dibaryon, albeit for 
non-physical values for the pion mass. NPLQCD's linear 



(quadratic) extrapolation to the physical pion mass gives 
B H = -0.2 ± 8.0 MeV (7.4 ± 6.2 MeV) @. Carames 
and Valcarce @ recently studied the H with a chiral 
constituent model constrained by AN, 'EN, EN and AA 
cross section data and find Bh values that are similar to 
the NPLQCD extrapolated values. 

These recent theoretical results motivate searches for 
the H with mass near the Mh = 2m A threshold. For 
masses below threshold, the H would predominantly de- 
cay via AS = +1 weak transitions to An, E°n or 
Ap7r~ final states. For masses above 2mA, but below 
m^o + m n (= 2mA + 23.1 MeV), the H would decay 
via strong interactions to AA 100% of the time. The 
E522 collaboration at KEK studied AA production in 
the 12 C(K~ , K + AAX) reaction and reported an intrigu- 
ing near-threshold enhancement but with limited statis- 
tics Q. The BNL-E836 collaboration searched for the 
AS* = +2 reaction 3 He(Ji ~, K + )Hn and established 
cross section limits spanning the range 50 MeV < Bh < 
380 MeV [IoT |. Searches for a bound H decaying to Apir~ 
reported negative results 



ll|,[l2j. Earlier searches, also 



with negative results, are listed in Ref. [13j. 

Decays of narrow T(nS) (n = 1,2,3) bottomonium 
(bb) resonances are particularly well suited for searches 
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for multiquark states with non-zero strangeness. The 
T(nS) states are flavor-5{7(3) singlets and primarily 
decay via the three-gluon annihilation process {e.g., 
(B(T(15) -> 3.93) = 81.7 ±0.7% [14]). The gluons mate- 
rialize into uu, dd and ss pairs in roughly equal numbers. 
The high density of quarks and antiquarks in the lim- 
ited final-state phase space is conducive to the produc- 
tion of multi-quark systems, as demonstrated by large 
branching fractions for inclusive antideuteron (D) pro- 
duction: B(T(15) -> DX) = (2.9 ± 0.3) x 10" 5 and 
B(T(25) DX) = (3.4 ±0.6) x 10" 5 [15]. An upper 
limit for the production of a six-quark S = —2 state in 
T(nS) decays that is substantially below that for the six- 
quark antideuteron would be strong evidence against its 
existence. 

Here we report results of a search for U-dibaryon 
production in the inclusive processes T(l,25) — > H X; 
H — » Apir~ and AA [16J. We use data samples con- 
taining 102 million T(1S) and 158 million T(25) de- 
cays collected with the Belle detector operating at the 
KEKB e+e~ collider [13] ■ The data were accumulated at 
center-of-mass system (cms) energies of y/s = 9.460 GeV 
and 10.023 GeV, which correspond to the T(1S) and 
T(2S) resonance peaks, respectively. Contributions from 
the e + e~ —> qq (q = u, d, s, andc) continuum pro- 
cess are inferred from a 63.7 fb _1 sample collected at 
y/s = 10.53 GeV and scaled by luminosity and 1/s. We 
assume equal T(1S) and T(2S) branching fractions: i.e., 
B(T{1S) -> HX) = B(T{2S) -> H X) = B(T(1,25) -> 
H X). 

Belle is a large-solid-angle magnetic spectrometer con- 
sisting of a silicon vertex detector, a cylindrical drift 
chamber (CDC), an array of aerogel threshold Cherenkov 
counters (ACC), a barrel-like arrangement of time-of- 
flight scintillation counters (TOF), and an electromag- 
netic calorimeter (ECL) comprised of CsI(Tl) crystals 
located inside a superconducting solenoid coil that pro- 
vides a 1.5 T magnetic field. Measurements of dE/dx in 
the CDC, ACC light yields, TOF flight times and ECL 
energy deposits are combined to form particle identifi- 
cation (pid) likelihoods C(h) (h = e + ,n + , K + ovp) for 
charged tracks. The K(h\h') = C(h)/{C(h) + C(h')) ra- 
tios are used to make pid assignments. Belle is described 
in detail elsewhere [18j. 

Samples of simulated T(IS') and T(2S) Monte Carlo 
(MC) events, generated with PYTHIA [19| and simulated 
using GEANT3 [23], are used to study backgrounds and 
determine efficiencies. For signal MC for various H de- 
cay modes, we use PYTHIA with the S*°(1530) mass, 
width and decay-table entries replaced with hypothesized 
parameters for the H. For MC-based optimization of 
selection criteria, we optimize a figure of merit defined 
as FoM = n S i g / y/n s i g + ribkg, where n S i g («bkg) is the 
number of selected signal (background) events assuming 
B(T(nS) -> HX) = 3 x lO" 5 . 

For both investigated channels, event selection starts 



with the identification of a A candidate reconstructed via 
its pn~ decay using the A-momentum-dependent crite- 
ria based on proton pid, track vertex information, decay 
length, and M(pn-) described in Ref. [Hj]. The M{pTt~) 
distribution for selected candidates is well fitted by a 
Lorentzian function with a FWHM resolution for the A 
peak of 1.50 ± 0.01 MeV. For A candidates, we require 
AM A = \M(jm) - m A | < 3.0 MeV. 

For the H — > Apir~ search, the pn~ track selection 
requirements are optimized using FoMs determined by 
MC assuming th = t\. Both the p and n~ are re- 
quired to be well identified by the pid measurements: 
K(p\h + ) > 0.9 (h+ = 7T+ or K+) n{*-\e-) > 0.9 
and K(n-\K-) > 0.6 0. We require that the p 
and 7T~ tracks and the A trajectory satisfy a fit to a 
common vertex with X\ p7r - < 50. In addition we re- 
quire CTx P K- > 0.0, where cr = I ■ phMh /\ph\ 2 and 
£ is the displacement between the run-dependent aver- 
age interaction point (IP) and the fitted vertex posi- 
tion. In some cases, the tracking algorithm finds two 
reconstructed tracks with nearly the same parameters 
from CDC hits produced by a single particle. Con- 
tamination from this source is removed by the require- 
ments M(pip 2 ) > 1878 MeV, M(?rf tt^) > 280 MeV and 
N h it s (pi)+N hits (p 2 ) > 50, where H -> Ap 2 ^2: A ~> Pi^i 
and A?hits (Pi) ^ s * ne num ber of CDC hits used to recon- 
struct the i th proton. In the Apir~ mode, there is a large 
background from A and p production via secondary in- 
teractions in the material of the beampipe and inner de- 
tector. This is removed by requiring \ph\ > 0.5 GeV for 
both h = A and h — p; this requirement is not applied 
to the Apir+ channel. In 6.3% (5.2%) of the data (MC) 
events, there are two or more entries that have one or 
more tracks in common. In these cases, the combina- 
tion with the smallest XA p7r - value is selected. For signal 
MC events, this chooses the correct combination 93.4% 
of the time. The A — > pi7rf candidate is subjected to a 
kinematic fit that constrains M(pi7rj~) to m^. The final 
selection efficiencies are determined from MC by aver- 
aging T(15) & T(25) signal MC to be ei = 7.7% for 
H -> Apn- and ei = 8.8% for H -> ApTT+ . 

The resulting continuum-subtracted Af(Ap7r~) 
(M(Apir + )) distribution for the combined T(15) and 
T(2S) samples, shown in the top (bottom) panel of 
Fig. [TJ has no evident H — > Apir~ [H — > Apir + ) 
signal. The curve in the figure is the result of a fit 
using an ARGUS-like threshold function to model the 
background [23j; fit residuals are also shown. 

For the second A (A 2 ) in the H ->■ AiA 2 (Aj -> PiK~) 
channel, in addition to the criteria used for Ai selection, 
FoMs based on MC events are used to optimize the addi- 
tional requirements Xa 1 a 2 < 200 from a AiA 2 vertex and 
IP constrained fit, and ct\ 2 > —0.5 cm. Entries in which 
two of the selected tracks originate from a single particle 
are removed by the requirements M(pip 2 ) > 1878 MeV, 
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FIG. 1: Top: The continuum-subtracted M(Apn~) dis- 
tribution (upper) and fit residuals (lower) for the combined 
T(1S) and T(2S) data samples. The curve shows the results 
of the background-only fit described in the text. Bottom: 
The corresponding A/(Ap7r + ) distributions. 



M(7rf ttJ) > 288 MeV and JV ut .(pi) + JV hi t s (p2) > 60. In 
3.2% (2.8%) of the data (MC) events, two or more en- 
tries have one or more tracks in common. In these cases, 
we choose the track combination that has the smallest 
XaiA 2 value. For signal MC events, this selects the cor- 
rect combination 95.4% of the time. The A A candidates 
are subjected to a kinematic fit that constrains both p?r~ 
masses to m\. The MC-determined selection efficiencies, 
obtained by averaging T(15) & T(25) signal MC results, 
are e 2 = 10.9% for H -> AA and e 2 = 10.1% for H ->■ AA. 

The difference between the AA and AA signal yields 
in the region M(AA) < 2.38 GeV, determined from two- 
dimensional fits to scatter plots of M(pi7ri) vs. M(p2^2) 
with the A mass requirements relaxed, is larger than the 
difference in the MC-determined AA and AA acceptances. 
This is attributed to deficiencies in the simulation of low- 
energy A and p inelastic interactions in the material of the 
inner detector. To account for this, a correction factor of 
0.83 ± 0.13 is applied to the H -> AA and H -> Kpir + 
efficiencies. The error on this factor is included in the 
systematic error. 

The continuum-subtracted M(AA) (M(AA)) distribu- 
tion for events that satify all of the selection requirements 
is shown in the top (bottom) panel of Fig. [5J where there 



FIG. 2: Top: The continuum-subtracted M(AA) distribu- 
tion (upper) and fit residuals (lower) for the combined T(15") 
and T(2S) data samples with the background-only fit super- 
imposed. Bottom: The corresponding M (AA) distributions. 



is no sign of a near-threshold enhancement similar to that 
reported by the E522 collaboration [§] nor any other ev- 
ident signal for H — > A A (H — > AA). The curve is the 
result of a background-only fit using the functional form 
described above; fit residuals are also shown. 

For each channel, we do a sequence of binned fits to 
the invariant mass distributions in Figs. [T] and [2] using 
a signal function to represent H —> /j (fi — Apir~ & 
f 2 = AA) and an ARGUS function to represent the back- 
ground. In the fits, the signal peak position is confined 
to a 4 MeV window that is scanned in 4 MeV steps across 
the ranges (toa + m p + m^-) < M(Apn~) < 2m^ and 
2m A < M (AA) < 2m A + 28 MeV. For the Ap7r~ (Ap7r+) 
mode, the signal function is a Gaussian whose resolution 
width is fixed at its MC-determined value scaled by a 
factor / = 0.85(1.12) that is determined from a compar- 
ison of data and MC fits to inclusive EE - — > Att~ and 
S°(2470) — > S~7r + signals found in the same data sam- 
ples. For the AA mode, the signal function is a Lorentzian 
with FWHM fixed at either T = or 10 MeV convolved 
with a Gaussian. Since the /, and fi acceptances are dif- 
ferent, we fit the particle and antiparticle distributions 
separately. 

None of the fits exhibit a positive signal with greater 
than 3cr significance. The fit results are translated into 
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90% CL upper limits on the signal yield, N^^Mh) and 
N^ l (Mh), by convolving the fit likelihood distribution 
with a Gaussian whose width equals the systematic error 
(discussed below) and then determining the yield below 
which 90% of the area above Nj, = is contained. These 
values are used to determine upper limits on the inclusive 
product branching fractions via the relation 

B(T(1,2S) -> H X) . B(H -> ft < 

1 N^(M H ) 

where N T = (260±6) x 10 6 is the total number of T(15) 
plus T(25) events in the data sample [2J| and B\^ p7r - = 
0.639 ± 0.005 Q. 

Sources of systematic errors and their contributions 
are listed in Table HI The tracking, pid and A reconstruc- 
tion uncertainties are common to other Belle analyses 
and are determined from data-MC comparisons of vari- 
ous control samples. For the channel-specific vertex re- 
quirements, we use data-MC differences found in high- 
statistics samples of inclusive T(l,2£) — > Apir + and AA 
events with M (Aprr) < 2.28 GeV (M (AA) < 2.38 GcV) 
selected with the same vertex criteria. The continuum 
subtraction systematic error contribution is determined 
from the errors in the relative on- and off-resonance lumi- 
nosity measurements. Systematic errors associated with 
the MC-determined acceptance and minimum momen- 
tum requirement are determined by varying parameters 
used in the PYTHIA generator and GEANT simulation 
programs. The systematic errors associated with the sig- 
nal fitting are determined from changes induced by vari- 
ations in the binning and fitting ranges in fits to an in- 
clusive S°(2470) — > S~7r + signal seen in the same data 
sample. Quadratic sums of the individual contributions 
are taken as the total systematic errors. 

TABLE I: Systematic error sources (in percent). When the 
H and H values differ, the H values are given in parentheses. 



Source 


H — » Apn 


H — > AA 


Nr(is) + Nr(2.s) 


2.3 


2.3 


tracking 


3.6 


3.6 


particle id 


7.2 


4.3 


A reconstruction 


3.2 (5.3) 


12.6 (9.6) 


Vertex requirements 


3.9 


3.5 


Signal efficiency 


2.0 (15.7) 


1.9 (15.8) 


Continuum subtraction 


1.4 


1.4 


B(A -> pn~) 


0.8 


1.6 


Fitting 


2.0 


2.0 


Resolution 


2.6 


2.6 


Quadrature sum 


10.2 (19.1) 


14.7 (19.8) 



For the final limits, we use the branching fraction value 
that contains <90% of the above-zero area of the product 
of the H and H likelihood functions. Figure [3] shows the 



resulting Mh — 2?7JA-dependent upper limits for the Ap7r~ 
and AA (for r — 0) modes. The upper limit values, listed 
in Table HI1 are all more than an order of magnitude lower 
than the average of measured values of B(T(1,25) — > 
D X), shown in Fig. [3] as a horizontal dotted line. 

The H — > Apir~ limits quoted in Table [TT| and shown 
in Fig.[3]are determined for an H lifetime tjj = 0.263 ns, 
i.e., the A lifetime. The acceptance decreases and, there- 
fore, the limits increase, with increasing lifetime: for 
th = 5ta, the acceptance is a factor of two lower and the 
limits are correspondingly twice as high. Conversely, for 
shorter lifetimes, the acceptance increases: for r = 0.5ta, 
the acceptance is higher and the limits are more stringent 
by 12 ± 2%. 
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FIG. 3: Upper limits (at 90% CL) for B(T(1, 25) -> H X) ■ 
B(H —> fi) for a narrow (F = 0) ff-dibaryon vs. Mh — 2m a- 
The vertical dotted line indicates the Mh = 2mA threshold. 
The limits below (above) the 2ttia threshold are for fi = 
Apir~ (/2 = AA). The horizontal dotted line indicates the 
average PDG value for S(T(1, 25) D X). 



TABLE II: 90% CL upper limits (xlO -7 ) on the product 
branching fraction B(T(1, 25) H X) ■ B(H -> fi), /i = 
Apir~; 6 Mi = 2mA - M H and f 2 = AA; 8M 2 = M H - 2toa. 



SMi (MeV) 


2 


6 


10 


14 


18 


22 


26 


30 


34 


fi = Apn~ 


15. 


9.7 


7.1 


6.3 


1.5 


5.2 


1.7 


4.6 


0.8 


f 2 =AA 




















r = o 


6.0 


9.6 


2.2 


11. 


14. 


9.2 


2.5 






r = 10 MeV 


16. 


17. 


15. 


37. 


44. 


42. 


33. 







The results reported here are some of the most strin- 
gent constraints to date on the existence of an H- 
dibaryon with mass near the 2mA threshold [25j]. These 
upper limits are between one and two orders of mag- 
nitude below the average of the PDG value for inclu- 
sive T(15) and T(2S) decays to antideuterons. Since 
T — > hadrons decays produce final states that are flavor- 
SU(3) symmetric, this suggests that if an ii-dibaryon 
exists in this mass range, it must have very different dy- 
namical properties than the deuteron, or, in the case of 
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Mh < 2mA, a strongly suppressed H — > Ap7r~ decay 
mode. 
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